Particle-hole bound states of dipolar molecules in optical lattice 



(N 

o 

(N 
D 

in 



CZ2 
I 

^— > ■ 

a- 

=5 



i 

O 

o 



cn 



X 



Yi-Cai Zhang 1 , Han-Ting Wang 1 , Shun-Qing Shen 2 , and Wu-Ming Liu 1 
1 Beijing National Laboratory for Condensed Matter Physics, 
Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China 
2 Department of Physics and Centre of Theoretical and Computational Physics, 
The University of Hong Kong, Pokfulam Road, Hong Kong, China 
(Dated: September 6, 2012) 

We investigate the particle-hole pair excitations of dipolar molecules in optical lattice, which can 
be described with an extended Bose-Hubbard model. For strong enough dipole-dipole interaction, 
the particle-hole pair excitations can form bound states in one and two dimensions. With decreasing 
dipole-dipole interaction, the energies of the bound states increase and merge into the particle-hole 
continuous spectrum gradually. The existence regions, the energy spectra and the wave functions 
of the bound states are carefully studied and the symmetries of the bound states are analyzed with 
group theory. For a given dipole-dipole interaction, the number of bound states varies in momentum 
space and a number distribution of the bound states is illustrated. We also discuss how to observe 
these bound states in future experiments. 
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I. INTRODUCTION 

Cold and ultracold molecules are attracting more and 
more attention due to their broad applications in the 
fields of high-precision measurement, quantum chemistry, 
quantum information and many-body physics |lj . Recent 
years the experimental techniques for cold and ultracold 
molecules have been developed greatly. The homonuclear 
molecular Bose-Einstein condensation (BEC) was real- 
ized experimentally 0, Q and the quantum state with 
exactly one molecule at each site of an optical lattice was 
created with a Feshbach resonance and STIRAP (stimu- 
lated Raman adiabatic passage) techniques [J, |5|. Mean- 
while the ultracold heteronuclear molecules were also pro- 
duced 

The heteronuclear molecules, such as SrO, RbCs, or 
NaCs, prepared in their electronic and vibrational ground 
states, have considerable permanent electric dipole mo- 
ment. The dipole-dipole interaction between molecules 
can be generated by the external applied electric field 
Moreover, the magnitude of the dipole-dipole inter- 
action can be controlled by the strength of the electric 

field pm. The tunable long range dipole-dipole inter- 
action may significantly modify the ground state and col- 
lective excitations of trapped condensates. For example, 
in trapped dipolar gases, it was shown theoretically |l3[ 
that the mean-field inter-particle interaction and, hence, 
the stability diagram are governed by the trapping ge- 
ometry. With increasing dipolar interaction, the ground 
state of rotating atomic Bose gases undergoes a series of 
transitions between vortex lattices of different symme- 
tries: triangular, square, "stripe", and "bubble" phases 
In rapidly rotating Fermion gas, the dipole-dipole 
interaction may even result in the fractional quantum 
Hall- like states [l5jj. 

The particles with long-range interactions in optical 
lattice can be described with extended Hubbard model 
[TB - [l8j . Compared with the regular Bose-Hubbard 



model with on-site interaction, the extended Bose- 
Hubbard model has richer ground state phases, such as 
Mott insulator, particle density wave, superfluidity or 
supersolid phase [19l426| . On the other hand, the gap- 
ful particle and hole excitations in the insulating phase 
[27h29I | may bind together and form bound states due 
to the long-range interactions. Although the excitons 
(holon-doublon pairs) in one-dimensional fermionic Hub- 
bard model have been extensively studied [3(|, the ex- 
citons in higher dimensions, and especially, in bosonic 
Hubbard models are less studied. 

In this paper, we study the particle-hole pair excita- 
tions of dipolar bosonic molecules in optical lattice, espe- 
cially, the possible bound states due to the dipole-dipole 
interaction. The paper is organized as follows. In Sec. 
II we introduce the extended Bose-Hubbard model to 
describe the polarized bosonic molecules in optical lat- 
tice, and derive the eigen equations to describe the sin- 
gle particle-hole pair excitation. In Sec. Ill, the ex- 
istence regions, the energies and the wave functions of 
the particle-hole bound states are calculated in one and 
two dimensions. The symmetries of the bound states are 
analyzed and possible experimental observation of these 
bound states is also discussed. A summary is presented 
in Sec. IV. 



II. THE PARTICLE-HOLE BOUND STATES IN 
THE EXTENDED BOSE-HUBBARD MODEL 



Considering only the nearest-neighbor interaction V to 
simulate the effect of dipole-dipole interaction, we write 
the extended Bose-Hubbard model as, 
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FIG. 1. Sketches of optical lattice and interactions of ultra- 
cold heteronuclear molecules, (a) The landscape of optical 
lattice potential in the x-y plane. The square optical lattice 
potential V„(x,y) = V op t[sin 2 (k x x) + sin 2 (k y y)] can be cre- 
ated by counter propagating far detuned laser beams, where 
V op t is the depth of lattice and k x r y ) is the laser wave vector 
along x(y) direction, (b) The interactions between dipolar 
molecules in the optical lattice (see equation(l)). The ar- 
rows on the molecules indicate the polarization of the electric 
dipole moments. 
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where np — b^b? is the number operator at site r with 
bp(bt) the annihilation (creation) operator of particle, a 
denotes the nearest neighbor vectors of site r. The on- 
site interaction U, the nearest neighbor interaction V and 
the hopping t are expressed as 



U = g I dx\W(x)\\ 

V = J dxdtf \W(x , -r)\ 2 V dd {x-x")\W(x* - (f+a))\ 2 , 

r h 2 
t = - / dxW*(x-r)(-— V 2 + V Q (x))W(x - (f+a)), 
J 2m 

(2) 

where W(f) is Wannier function corresponding to the 

lowest energy band and Vdd(%) = iff 1- fg|° 3 8 1S tue 
dipole-dipole interaction of two particles separated with 
a distance \x\. For the electric dipole-dipole interaction, 
Cdd = fjj- with d the electric dipole moment, eo is vacuum 
permittivity. We assume that the polarization of dipole 
moment is along the z axial direction. In the following 
studies, we consider particles confined in one-dimensional 



chain or two-dimensional square lattice and treat hopping 
term H\ as a perturbation. 

In the atomic limit (t = 0), n? is a good quantum 
number and the cigenstates of Hq can be written as direct 
product of number states. When the filling factor is one, 
V > and zV < U, the Mott insulating ground state can 
be written as 1 11111111...) with the ground state energy 
Eq = t;NzV, where z is the coordination number and N 
the number of lattice sites. The excited states contain 
one or more particle-hole pairs. A single particle-hole 
pair state with a hole at r and a particle at r' can be 
represented as \r,r') — |110 r 111112 r /ll...). Taking into 
account the hopping term t, the particle and hole will 
move in the lattice. Similar to the two-magnon states 
in ferromagnetic system |31| , the single particle- hole pair 
state can be written as a linear combination of \r, r'): 



\*l>) = E 



r, rf) 



(3) 



where (/)?,?' will be determined by solving the approxi- 
mate eigen equation in single particle-hole pair subspace 



(f,ff\H\ip) = E{r,r/\i/j). 



(4) 



Calculating H\f, r 1 ) and considering the boundary con- 
dition (j>ff ff — (a particle and a hole do not share the 
same lattice site), we get 



+ V 5 P ^ s [V^ r -, + {-m ?+3 p + {-2t)4> ? -»_ J 



— <5r'-r,0 ^^(^0f+ff,f / + 2t(j>ff + ff), (5) 
B 

where 8pt—? a is the Kronecker delta function and e = 
E — Eq — U . The particle-hole excitation energy is u — 
E-E = e + U. 

Owing to the translational invariance, it is convenient 
to apply a transformation — '^ elK R< l ) {p)i where 
R = and p = r — r denote the center-of-mass and 
the relative coordinates respectively and K the total mo- 
mentum of the particle and hole. At fixed K, the eigen 
equation (5) becomes 



- Y,[te i& ^<j>{p- a) + 2fe'^^(p + «?)] = 

a 

e^p) + hAV<t>(p) + (~te^ - 2te^ )4>{p- a)] 

- <V,o5>e-^ <f>(a) + 2te iS ^4>{5)]. (6) 



Considering the finite range character of the nearest 
neighbor interaction V, we utilize the Green's function 
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FIG. 2. Energy spectra of the bound states of particle- 
hole pair in one dimension for different nearest-neighbor in- 
teractions V — 2t (black circles), V — 3t (red triangles), 
V — 5t (black diamonds) and V = 8t (red squares) respec- 
tively (K = K x ). As comparison, the continuous spectrum 
of a particle-hole pair excitation is also shown (green shaded 
region). 



approach to solve the eigen equation. To this end, we 
introduce an effective Hamiltonian 



where 



Hph,a — — t ^ ( 



,iK-ff/2 



\p)(p-a\-2te 



iK-ff/2 



(7) 



\P)(P + Z\, 



& 

a 

to describe the motion of a particle around a hole, where 
\p) is basis set in relative coordinate spaces. The in- 
dependent H p h.Q describes the kinetic energy of the par- 
ticle and Vph denotes the interaction between particle 
and hole. With \<j>) = ^2^^>(p)\p), it is easy to reproduce 
Eq. ([6]) from the eigen equation H e ff\<f>) = e\4>). After a 
Fourier transformation, the kinetic energy is obtained as 



i{K/2-q)-3 



, (8) 



FIG. 3. Probability density of a bound state in real space 
in the case of V = 8t and K = 0. The probability density 
decreases exponentially with increase of distance. 



Introducing retarded Green functions Go = 
linr ?? ^o+ l/(e — H p h,o + irj) and G = lim^_j.o+ — 
H e ff +if]), we could calculate G through the Lippmann- 
Schwinger equation G = Go + GoV p hG. In the real 
space, the Green's function G has a general form 



(p'\G\p) = limV 



<M/?')</4(p) 

e - E m + irj 



G{p,p,e), 



where the contributions of the continuous spectra are ne- 
glected and E m 's are the eigenvalues of the bound states 
with (j> m 's the corresponding eigenfunctions [32| . We can 
determine the eigenvalues and eigenfunctions by analyz- 
ing the poles and residues of the obtained Green's func- 
tion. 

From Eq. (J7]), we have 

.-}.- _ ,. 1 d d qe^<P'-^ 

{p'\Go\p}= am — — ^ / -, 

, ; ^o+ (27r) d J_ n £ - £q [K, q) + ir) 



and 



(p'\G\p) = (p'\G \p) 

+ J2 {p'\Go\p-i){p-i\Vph\P2){p2\G\p). (9) 



where q is the relative momentum of particle and hole. 
When a particle is adjacent to a hole, the energy offset 
— V in Vph may result in the formation of particle-hole 
bound states. 



For specific pi — or a, Eq. (J9]) can be reduced to 
a set of simultaneous linear equations. Green's func- 
tions (p'\G\p) with pi = or a can thus be obtained 
exactly. The residues of (Q\G\p) is always vanishing for 
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FIG. 4. The mean size d of bound states as a function of V at 
K=0 in one dimension. When the nearest-neighbor interac- 
tion V approaches the critical value V — 3t, the mean size of 
the bound state becomes divergent and the bound state ener- 
gies merge into the continuous spectrum (see the red triangles 
in Fig. 2), indicating the disintegration of the bound states. 



non-vanishing bound states energies e in our calcula- 
tions, which is consistent with the boundary condition 
(j)(0) = 0. We present our results in the next section. 



III. BOUND STATES IN ONE AND TWO 
DIMENSIONS 

A. Bound states in one dimension 
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In the Mott insulating phase, the excitation energy of 
a single particle/hole was calculated as 

Sq 1 

u p . h = ±[-— + U(n Q - -) + 2dVn Q - fi] 



+ [{^? + eoU{n Q + 1 -) + ^ (10) 

with the dynamical Gutzwiller approach 29]. Up to the 
first order of jj, the single particle-hole pair excitation 
energy is 



lo = uj p (k p ) + Lu h (k h ) = U + e (K, q) 



(11) 



with the total momentum K = k p + kh and relative mo- 
mentum q= (kh—k p )/2. For the extended Bose-Hubbard 
model, however, a particle-hole pair may form a bound 
state due to the interaction (See Eq. ( ffj l). In the follow- 
ing calculations, we take the lattice constant as 1 and use 
t > as the unit of energy. Conditions oit/U <C 1 and 
zV < U are assumed to make sure that the ground state 
of the system is the deep Mott insulating phase. 

In one dimensional case, So(K,q) = — 2tg cos(g + 

6 X ) with t & = t^9cos 2 (^) + sin 2 (^) and 6 & = 



FIG. 5. The wave functions of bound states in two dimen- 
sional case at K = (0,0) and V = 12t (non-normalized), s-, 
p x - and d a .2_ !/ 2-wave symmetries can be observed respectively. 



arctan(jjtan(^f-)) . The top and bottom boundaries of 
particle-hole continuous spectrum are 



(12) 



(see Fig. 2). Compared with particle-hole pairs in fermion 
Hubbard model with half filling [34j, the minimum of 
band width of the continuous spectrum is non-zero at 
K = ±7r. This is because that there are no particle-hole 
symmetries in Bose-Hubbard model. 

For a specific nearest-neighbor interaction V and total 
momentum K, we search for bound state solutions out- 
side the continuum. In the case of e < —2t±, the free 
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FIG. 6. Variations of the energies of different bound states in 
two dimension at K = (0, 0). The energies of different bound 
states increase with decreasing V and cross the particle-hole 
continuum (— \2t at K — (0,0)) at different V. 



Green function is calculated as [33J 

1 r d d qe^<P'-^ 



(p'|G |p)= lim 



(13) 



whe 



re a = — 



= , 7 = e %Bss and (3 = e K with K = 



arccosh^T^:). 

With p' = 0, ±1, Eq. © is reduced to 3 linear equa- 
tions. The Green functions of (— l|G|p), (0|G|p) and 
(l|G|p) can be exactly obtained as 



Q{-p){P/i)- {p+1) 

-Utx-V/3) 



(-l\G\p) 
<0|G|p>=0, 
(HGIp)^)^- 1 



i(tz-Vf3) 



p^0 



(14) 



where 9{p) is Heaviside step function. From these Green 
functions, we obtain two degenerate bound states corre- 
sponding to the particle on the left (right) of the hole 
respectively. The bound state energy is 



t 2 . 
V 



(15) 



Accordingly, j3 is simplified as /3 = < 1 and the con- 
dition for the existence of the bound state is obtained 
as 



(16) 




FIG. 7. Dispersion relations of bound states along the line of 
K x = K y — K for V = 6t and V = 12t. As comparison, the 
particle- hole continuum is also shown (green shaded region). 



When < V < t, the interaction between particle and 
hole is too weak to bind them together and no bound 
state is found. When t < V < 3i, two degenerate bound 

states exist in the regions of 2arccos\ 



v 2 -t 2 



< KJ < IT. 



The region get larger with the increase of V. When V > 
3t, bound states may be found in the whole first Brillouin 
zone. In Fig. 2, we show the spectra of the bound states 
for the interaction V = 2t, 3t, 5t and 8t, respectively. 
The energies of the bound states decrease with increasing 
V. At V = 2t, the calculated existence intervals are 

0.587T < lifa-l < 7T. 

From the residues of the Green function, we can also 
extract the bound state wave functions. For example, 
the bound state with a particle on the right of a hole is 
written as 



MP) = cfpr-Ve-^'-V p >l 

4>r(p) = o P < l, 



(17) 



t& is the normalized constant. In 



where C = \l 1 
Fig. 3, we show the probability distribution of the wave 
function for K = and V = 8t. 

The mean size of the bound state is calculated as 
d = ((/>r(p)\p\Mp)) = yT~pr ■ In Fig.4, we show the 
mean size of the bound state as a function of the nearest- 
neighbor interaction V at K x — 0. The stronger is the 
interaction, the smaller is the size and the closer do the 
particle and hole bind. 
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B. Bound states in two dimension 

Before presenting the numerical results, we discuss the 
symmetries of the effective Hamiltonian H e f f in Eq. ([7]) 
in details. After a gauge transformation 



\p)' = \x,y)' = e 



-i{x9%-\-y9i 



)\x,y) = e-^+^lp) 



(18) 

with 9 X = arctan(-^tan(^f L )) and 9y=arctan(^tan(-^ L )) , 
we could remove the phase factors in the Hamiltonian 
and get 

H 'eff = H 'ph,0 + V ph> ( 19 ) 
H 'ph,0 = ~*i \P)'(P~ °*\' - h W{P- S v\'i 

Vph = - v Yl + ** E i^>'<°i' + ^ 

where S x and er y denote the nearest neigh- 
bors along the x and y directions, respec- 
tively, and t£ = tJ§cos 2 (^f) + sin 2 (^f) and 



t ,) 



t\J Qcos 2 (^f-) + sin 2 (^ L ) are effective hoppings 
along the x and y directions. 

When K x — ±K y , t x = t y , H' e ff has symmetries of D4 
group. According to the irreducible representations of 
this group, the bound states can be classified and labeled 
with s, p x ( y ) and d x 2_ y 2 wave, respectively. Among them, 
the s wave belongs to an identical representation A\ of 
D4 group, the degenerate p x and p y waves form a two di- 
mensional irreducible representation E of D4 group, and 
the d x 2_ y 2 wave belongs to an irreducible representation 
B\ of D4 group. 

When K x 7^ ±K y , t x 7^ t y , the symmetry reduces to 
D2, a subgroup of D4. For simplicity, we still label the 
four bound states with s, p x (y) and d x 2_ y 2. Differently, 
here all the s and d wave belong to identical represen- 
tations A\ of D2 group, while p x (Py) wave belongs to 
an irreducible representation B^{B2) of D% group. The 
degeneracy of p x and p y waves is lifted. 

With p' = (0,0), (±1,±1), Eq. (9) reduces to five lin- 
ear equations in two dimensional case. The free Green 
function Go(p, (7) can be expressed with elliptic integrals 
[HI, [35[. Although no brief solutions of < p\G\p' > could 
be found, we may factorize the particle-hole bound state 
equations and analyze the existence conditions for the 
bound states along the symmetric lines of K x = ±K y = 
K [35] . Considering the asymptotic behaviours of elliptic 
integrals, we get the thresholds of the interaction V as 
follows: V cr ,s(K) — t x for s wave, V cr ,p{K) = 2 _ 4 ta for 
p x ( y ) wave and V cr ^d{K) — 8 _2 7r ^a f° r d x 2_ y 2 wave. From 
these thresholds, the existence region of every bound 
state is obtained respectively. For s wave, we have 



<K 



2arccos 



V 2 - 1 2 



When V > 3t, the s wave bound state may be found in 
the whole Brillouin zone. The existence region of p wave 
bound states is 



2arccosJ <\K\<«,Ve (^, ^t]. 



When V > -^si, the existence region extends to the 
whole Brillouin zone. For d wave, existence region is 
expressed as 



2arccos\ 



8t 2 



< \K\ < tt,V e ( 



3tt 



4 - 7T 4 - 7T 



8/ 2 



< \K\ <ir,V e (i,3t] 



When V > -r^-t, the d wave existence interval is the 
whole Brillouin zone. 

Now we present the numerical results in two dimen- 
sions. At K = (0,0), four bound states could be found 
when V > V cr ^{0) = 10.97t. The wave functions with 
s— , p x — and d x 2_ y 2— wave symmetry for V = I2t are 
shown in Fig. 5. With the decrease of V, the high- 
est d x 2_ y2 -wave, the degenerate p^^-waves and the low- 
est s-wave merge into the continuous spectrum one by 
one. Finally, all the bound states disappear when V < 
Vcr,s{0) = In Fig. 6, we illustrate the variations of 
the bound state energies with the changes of V. 

We then search for bound state solutions along the 
line of K x = K y = K. As shown in Fig. 7, there are 
four bound states for all K £ [— tt, tt] at V — I2t. At 

V = 6t, s-wave exists in the whole region, while p x ( y )- 
and d x 2_ y 2-wave states appear at 0.52-7T < \K\ < n and 
0.697T < I if I < tt respectively. Decreasing V further, 
we find that the regions of the bound states shrink, and 
all the bound states disappear when V < V cr ^ s (±7r) — 
It. Similar results are obtained along the line of K x = 

-Ky = K. 

Away from the lines of K x — ±K y , the degeneracy of 
p x and p y bound states is lifted, as mentioned before. At 
a given V , different bound states have different existence 
regions in K space. Consequently, the number of bound 
states vary in the (K x ,K y ) space. As an example, we 
show the number distribution of the bound states for 

V = 6t in Fig. 8. 



C. Discussion of observations of the bound states 

Inelastic light scattering directly measures the dynami- 
cal structure factor S(q, uj), the Fourier transformation of 
density correlations. Bragg spectroscopy has been pro- 
posed to detect quantum phases in optical lattice and 
successfully applied to measure the excitation spectra 
(phonons of BEC), the composition of the excitations and 
the Higgs-type amplitude mode in the superfluid conden- 
sate, as well as the particle-hole excitation energies in the 
Mott-insulator state [3614421 ] . When this technique is uti- 
lized in the study of ultra-cold polar molecules in the 
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FIG. 8. The number distribution of bound states in the first 
Brillouin zone at V = 6t. As shown, there are more bound 
stats at the corners than that at the center in the first Bril- 
louin zone. 



optical lattice, we may expect extra resonance peaks cor- 
responding to the particle-hole bound states lying outside 
the particle-hole continuum. 

Compared with the traditional solid state counter- 
parts, the deep Mott insulating state with t/U <C 1 
and zV < U can be realized by tuning the inter- 
action parameters t, U and V in ultra-cold dipolar 
molecules in the optical lattices. In deep optical lat- 
tice the hopping is approximately evaluated as t — 

(i/^EriVopt/Erf/^Xpi^Vopt/Er) 1 / 2 }, with V opt = 

18 ~ 20E r the optical lattice depth and E r the lattice re- 
coil energy 43] . Under the condition of — <C 1 , where a s 



is the s-wave scattering length and a± 



with uj± 



transverse trapping frequency in one dimension or trap- 
ping frequency of the z direction in two dimensions, the 

on-site interaction is estimated as U\d = ^J^fy^±a s /io 
and U 2D = ^flh with = (#-)*- and a the lattice 

mar( 2ir)li u v V opt ' ir 

constant [43l. |44||. 

Taking Bose molecule 23 Na 7 Li for example, we may 
estimate the typical values of t, U and V for observing 
the particle-hole bound state. The molecule prepared 
in the ground state has a permanent electric moment 
d = 0.58D [451] . With a lattice constant a ~ 0.5 [im, 



a s /a ~ 0.01 and the transversal tapping frequency cjj_ ~ 
2n x 10 4 Hz, we have t ~ InK , U\d ~ 25nK and U2D ~ 
50nK. The nearest-neighbor interaction V can be tuned 
as V ~ WnK with an applied electric field. 



IV. SUMMARY 

In summary, we have investigated bound states of 
particle-hole pair resulting from the dipole-dipole inter- 
action between polar molecules in the optical lattice. For 
a large enough dipole-dipole interaction, two degenerate 
bound states, which correspond to a particle on the left 
and the right of a hole, are shown to exist in one dimen- 
sion. While in two-dimensional case, four bound states, 
with s-, p x (y)- and d x 2_ y 2- symmetry respectively, are 
found along the lines of K x = ±K y . Away from the lines 
of K x — ±K y , the degeneracy between p x and p y waves is 
lifted. With decreasing the nearest-neighbor interaction 
V, the energies of the bound states increase and merge 
into the particle-hole continuum gradually. The wave 
functions, the dispersion relations and the existence re- 
gions of the bound states are studied in details. For a 
given nearest-neighbor V, the number of bound states 
is different in different regions of K space and a num- 
ber distribution of bound states is given for the nearest- 
neighbor V — 6i in two dimensions. The possible exper- 
imental observation of the particle-hole bound states is 
also discussed. 

Electron-hole bound state excitations (excitons) have 
been extensively studied for many years and very re- 
cently, the excitonic Bose-Einstein condensates have been 
realized experimentally in semiconductors [46l . |47| . We 
hope our study on the particle-hole bound states in 
bosonic systems would enrich our understanding of ele- 
mentary excitations in quantum many-body systems and 
stimulate more efforts on the bound state phenomena in 
various systems such as magnets, superconductors and 
atomic systems. 
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